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We calculate the superconducting phase diagram as a function of temperature and z-axis 
anisotropy in a model for tetragonal and cubic systems having strong antiferromagnetic fluctua- 
tions. The formal basis for our calculations is the fluctuation exchange approximation (FLEX) 
applied to the single-band Hubbard model near half-filling. For nearly cubic lattices, two supercon- 
ducting phase transitions are observed as a function of temperature with the low-temperature state 
having the time-reversal symmetry-breaking form, d^2_y2 ii id^^2_^2. With increasing tetragonal 
distortion the time-reversal-symmetry-breaking phase is suppressed giving way to only d^2_y2 or 
dsz^-r^ single-component phases. Based on these results, we propose that Celus is a candidate for 
exhibiting a time-reversal symmetry-breaking superconducting state. 



The discovery of unusual superconductivity in the 
heavy fermion systems focused much early effort on elec- 
tronic pairing mechanisms and the possibility that the 
resulting superconducting states have lower symmetry 
than the underlying crystalline lattice. The competi- 
tion between superconductivity and magnetic ordered 
states suggests that spin fluctuations are a prime can- 
didate to form the glue that binds electrons into Cooper 
pairs. Power law dependencies of thermodynamic prop- 
erties and phase diagrams that include multiple super- 
conducting states, for example the H — T phase diagram 
of UPts and the x — T phase diagram of Ui_2:Tha;Bei3 
point to the lower symmetry of the superconduct- 
ing order parameter. As a consequence of having a lower 
symmetry in comparison to that for the lattice, the order 
parameter of an unconventional superconductor can have 
more degrees of freedom than exhibited in a conventional 
superconductor leading to phase diagrams involving mul- 
tiple superconducting states, but microscopic models are 
needed to shed light on the connection between the pair- 
ing interaction and the pairing states that are produced. 

Recently, the interplay between crystal lattice sym- 
metry and spin-fluctuation induced pairing instabilities 
was explored by Monthoux and Lonzarich |3] and Arita, 
Kuroki and Aoki They find that spin fluctuation- 
induced pairing into the unconventional (1^2 _y2 state is 
most effective for producing large transition tempera- 
tures (Tc) for quasi two-dimensional lattices, a result 
that is consistent with the observation that the highest 
Tc values occur in the quasi two-dimensional cuprates 
and the more recent finding that Tc in the cubic heavy 
fermion system Celna 0| is about one order of mag- 
nitude less than is observed in a collection of related 
quasi two-dimensional compounds, such as the series 
Ce„TmIn3„_|_2m where T ~ Rh or Ir and n = 1 or 2 and 
771 = 1 5]. 

In this Letter we report on model calculations of the 
superconducting phase diagram as a function of tem- 



perature and tetragonal anisotropy for electrons paired 
via spin fluctuations in nearly antiferromagnetic systems. 
The phase diagram is obtained numerically through the 
proper generalization of the fluctuation exchange approx- 
imation for the single band Hubbard model to the super- 
conducting state, a generalization which is necessary to 
resolve the relative stability of nearly degenerate uncon- 
ventional pairing states. We find that for cubic lattices 
the stable superconducting state has the time-reversal 
symmetry breaking form d,j.2_y2 ±1(1^^2 _^2 . Small tetrag- 
onal distortions lift the degeneracy between 6,^2 _y2 and 
(i322„r2 pairing states leading to two superconducting 
phase transitions as a function of decreasing tempera- 
ture, but as the degeneracy is lifted further through larger 
tetragonal distortions the second transition is suppressed. 
On the basis of these results, we suggest that the low- 
temperature superconducting state in simple cubic sys- 
tems with strong antiferromagnetic correlations, such as 
Celua, are strong candidates for realizing multiple su- 
perconducting phases including a low-temperature phase 
with broken time-reversal symmetry. 

The microscopic basis for our calculations is the single- 
band Hubbard model, 

= - ^ {tij c\ „Cj,„ -f h.c.) + ^ XI cl.T^^T^i'^^d: 

(1) 

where tij represents the electron hopping amplitude be- 
tween sites i and j and U is the on-site interaction energy 
between up and down spin electrons. The values of tij 
reflect the underlying lattice structure. For tetragonal 
lattices, the hopping parameter for unit displacements 
in the x or y directions, t,j.y^ is distinct from the same 
for unit displacements along the z axis, t^- For simplic- 
ity we set the hopping integrals equal to zero for larger 
displacements in which case the non-interacting electron 
bandwidth is equal to Wq — 8txy + itz ■ 

An approximation scheme must be employed for cal- 
culations based on this model. We use the fluctuation 
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exchange approximation (FLEX) of Bickers, White, and 
Scalapino [Q], a numericaUy-based scheme that is con- 
serving in the sense described by Baym FLEX pro- 
vides a self-consistent description of both the quasiparti- 
cles and the magnetic-fluctuation induced pairing inter- 
action for a given on-site interaction strength, tempera- 
ture, and band fiUing. It is expected that results obtained 
with FLEX are quantitatively accurate for, at best, weak- 
to-intermediate couphng, i.e. U/Wq ^ 1 However, 
even within this range there are notable qualitative fail- 
ures such will be discussed shortly, violations of 
the Mermin-Wagner-Hohenberg theorem for c?-wave su- 
perconducting states. 

In contrast to other works, our formulation of FLEX 
includes the entire set of fluctuation diagrams for the 
electron self-energy in both the normal and supercon- 
ducting states, i.e. we include all possible combinations 
of particle-like, hole-like and anomalous Green's func- 
tions. At the formal level, this ensures that response 
functions and Green's functions are obtained in a consis- 
tent manner such that conservation laws derivable from 
symmetries of the Hamiltonian are obeyed. While Tc is 
determined by self-energy diagrams which contain only 
one anomalous Green's function, a more complex set of 
diagrams with three anomalous Green's functions con- 
tribute to fourth order terms in the Ginzburg-Landau 
expansion for the free energy in terms of the supercon- 
ducting order parameter and these terms are known to 
determine the relative stability of multicomponent pair- 
ing states 0. We note that the results we obtain for Tc 
in two dimensions are in agreement at the 10% level with 
those obtained earlier 0| suggesting that self-energy di- 
agrams that are omitted in most treatments of FLEX do 
not play a large role in determining Tc, at least in the 
two-dimensional limit. 

To access sufficiently large system sizes with modest 
computational resources, we combine this formulation of 
FLEX with the dynamical cluster approximation (DCA) 
[lll |. Essentially, in the DCA large lattice calculations 
are made feasible by approximating correlation effects 
via a smaller embedded cluster. When combined with 
the DCA, there are three numerical parameters in FLEX: 
the number of Matsubara frequency points used (m), the 
lattice size {Nl) and the DCA cluster size {Nc). We show 
the dependence of results for the pairing amplitude, nip 
(defined below), versus temperature on these numerical 
parameters in Figure^for a two-dimensional lattice. For 
the range of these parameters that we can access, the 
most significant variation in these curves is through the 
DCA cluster size (bottom graph in Figure QJ. Nonethe- 
less, we find that the error made in the result for Tc is only 
on the order 10% when using a 4^ DCA cluster. In what 
follows, we use m = 16384 frequency points, Nl = 32'^ 
lattices and Nc = 4^ DCA clusters. 

Although a variety of electronic and thermodynamic 
properties can be calculated with the FLEX, our main 



focus here is on results for the superconducting tran- 
sition temperature and the associated order parameter 
symmetry. We obtain these results by determining the 
spatial dependence of the pair wave function, "0(1"), for 
the stable superconducting state and the associated pair- 
ing amplitude, nip. These quantities are obtained from 
the self-consistent result for the anomalous Green's func- 
tion, F, via 

PitiT 0",r) = (cr=o,TCr,i) = V(r). (2) 

Other than requiring that tp{r) is even as a function of 
r {i.e. we restrict ourselves to singlet-pairing), no as- 
sumption is made on the symmetry of the pairing state. 
To permit any possible singlet-pairing state to emerge in 
these calculations, we initialize the self-consistent FLEX 
equations with a small, spatially random pairing field 
to induce a pairing amplitude in all possible symmetry 
channels. The small field is removed as the self-consistent 
procedure projects the wavefunction of the most stable 
pairing state, i^ir). 

We emphasize the mean-field nature of the FLEX 
phase diagram for d-wave superconductivity in the single- 
band Hubbard model. The mean-field Tc values essen- 
tially represent the temperature at which superconduct- 
ing order emerges locally. The true thermodynamic Tc 
is determined by phase fluctuations that are not present 
in FLEX. The neglect of these fluctuations is especially 
problematic in the two-dimensional limit where they are 
known to eliminate the possibility of a flnite tempera- 
ture phase transition. Nonetheless, it is well known that 
a relatively weak interplanar coupling can stabilize su- 
perconductivity in a quasi-two-dimensional systems and 
the mean-field phase diagram is revealing with respect 
to determining the conditions under which the tendency 
toward superconducting order is greatest. 

Our primary result, shown in Figure [3 is the calcu- 
lated superconducting phase diagram as a function of 
the scaled temperature, fc^T/Wo, and the ratio of in- 
terplanar to intraplanar hopping, t^/txy for fixed values 
of density {n — 0.85 electrons per site) and the ratio 
of the on-site interaction energy to bare electron band- 
width {U/Wq = 0.5). We focus on three features in Fig- 
ure El the low-temperature stability of the time-reversal 
symmetry breaking state dx2_y2 ± id^z'^_^2 for lattices 
with cubic and nearly cubic symmetry, two distinct su- 
perconducting transitions as a function of temperature 
for nearly cubic lattices, and, in agreement with previous 
results, maximal Tc values occurring in the quasi- two- 
dimensional limit. 

We first focus on the cubic limit, tz/txy — 1, for which 
Tc is at a local minimum. We find that, in agreement 
with Arita, Kuroki and Aoki that the stable super- 
conducting state belongs to the two-fold degenerate rep- 
resentation Fj^ of the cubic group which is described by 
basis functions of the form dx2_y2 and (i322_,.2. Sym- 
metry considerations allow for either single-component 
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or multicomponent states of the form dx2_y2, d^z2_^2, 
dx2-y2 id^z^-r^, or dx2_y2 iid^z^-r^. Weak-coupling- 
based arguments suggest that the most stable state is 
the one for which the superconducting gap is most com- 
plete on the Fermi surface as this will tend to maxi- 
mize the condensation energy. This argument favors the 
d^2_y2 ±id^z^-r^ pairing state as it has point nodes while 
the others have line nodes [l^ . 

FLEX incorporates feedback effects between quasipar- 
ticles, the order parameter and the pairing interaction 
and, thus, does not necessarily generate the pairing state 
expected from weak-coupling theory. Nonetheless, we 
indeed find that FLEX produces the dx2_y2 ± «d3j,2_^2 
pairing state in this case. This state breaks time re- 
versal symmetry leading to unusual phenomena such as 
as bulk magnetic effects associated with the supercon- 
ducting pairs U . Assuming that quasiparticles are well- 
defined in the superconducting state, such pairing is ex- 
pected to generate thermodynamic properties that reflect 
the point node structure of the gap function, such as hav- 
ing a low-temperature specific heat. 

For slight deviations from cubic symmetry, i.e. 
tz/txy = 1 i e, e ^ 1, increases due to the improved 
stability of either the dx2_y2 or d^z'^_r2 pairing state. 
The two pairing states are no longer degenerate for e 7^ 
and the initial transition from the superconducting tran- 
sition is into a single component state. However, because 
of the near degeneracy of these states, a second super- 
conducting transition occurs into the dj;2_y2 + id^^z'^-r^ 
state with decreasing temperature. The low-temperature 
state is nodeless because d^.z'^-r'^ basis functions are mem- 
bers of the symmetric group for tetragonal lattices and, 
consequently, are intrinsically mixed with contributions 
from basis functions with s-wave symmetry. However, s- 
wave terms in the pair wavefunction are relatively small. 
For example, the s-like terms have amplitudes which are 
about 4% of those for d^^z'^-r^ for tz/t„j = 1.1. Con- 
sequently, the gap is correspondingly small for a set of 
points on the Fermi surface that correspond to the point 
nodes that are found for cubic lattices. 

As is expected, we find that for lattices favoring in- 
plane conduction, i.e. < tz/txy < 1, the dx2-y2 pair- 
ing state is most stable, while the d^z^-r^ (+ s-wave) 
pairing state is most stable when interplanar motion is 
enhanced, i.e. tz/t^y > 1. ft is apparent from Figure [21 
that the largest Tc values are obtained in the quasi-two- 
dimensional limit, i.e. tz/trn—^ 0, in agreement with the 
results referred to earlier [2, y| and in line with the trend 
observed in the compounds Celua 0| and Ce„TmIn3„4.2m 
where T — Rh or Ir and n = 1 or 2 and m = 1 [5j. 
There also is a region of enhanced Tc for tz/t^y > 1 up 
to approximately tz/txy 2 at which point Tc becomes 
vanishingly small. We note that tz/txy = 2 corresponds 
to the point at which the interplanar bandwidth equals 
the bandwidth corresponding to in-plane motion. 

It is interesting to consider these results in light of 



experimental data for superconductors with cubic sym- 
metry. Data for the compound Pr04Sbi2 is consistent 
with the superconducting phases observed in this cal- 
culation [1^ . For example, specific heat measurements 
show evidence of a double superconducting transition 0| 
and muon-spin relaxation data is consistent with a low- 
temperature time- reversal symmetry breaking state 'l5^ . 
However, the strong collective mode that forms the glue 
between paired electrons is due, most likely, to electric 
quadrupolar rather than magnetic degrees of freedom 
[l6| . Nonetheless, it is interesting to note the possibil- 
ity that two different pairing mechanisms may tend to 
generate the same time-reversal symmetry breaking state 
in cubic symmetry and, thus, the phenomenolo gv d evel- 
oped for the superconducting state in Pr04Sbi2 [13 may 
apply also for magnetically-paired superconductors. 

The alloy series Ui-ajTh^jBeis displays a complex su- 
perconducting phase diagram as a function of thorium 
concentration, x. For 0.2 ^ a; ^ 0.4 values, a double 
superconducting transition is observed in the electron 
specific heat and muon spin resonance data points 
to the appearance of an internal magnetic field below 
the second transition. Sigrist and Rice interpreted the 
magnetic anomaly in terms of a multicomponent time- 
reversal symmetry breaking superconducting state |l8l |. 
However, Kromer et al. [19| use specific heat and lat- 
tice expansion data to argue for the appearance of spin 
density wave below Tc. The interpretation of the phase 
diagram for these compounds is still being debated [ioj . 
We simply note that the scenario described by Sigrist and 
Rice is consistent with pairing states that are observed in 
the model that we consider here and should their inter- 
pretation be incorrect it would point to the importance of 
processes that are neglected in spin- fluctuation models. 

The cubic superconductor Celna is, perhaps, the most 
likely candidate system for connecting with these model 
calculations. On account of the vicinity of the Neel state, 
there is a strong likelihood that electron pairing is related 
to the strong antiferromagnetic correlations in this sys- 
tem. To the best of our knowledge, no evidence has been 
presented either for a double superconducting transition 
as a function of temperature nor for a low-temperature 
time reversal symmetry breaking state in this compound. 
In light of weak coupling arguments jl2|, the FLEX re- 
sults presented here and the anomalies observed in other 
unconventional cubic superconductors, the presence or 
absence of time-reversal symmetry breaking in this com- 
pound will be revealing with respect to the applicability 
of simple models, such as the one considered here, for 
providing a minimal microscopic basis for understanding 
superconductivity in Celns and related compounds. 

In summary, we have performed fluctuation exchange 
approximation calculations for the Hubbard model near 
half-filling to model the phase diagram of tetragonal and 
cubic superconducting systems whose pairing is mediated 
by antiferromagnetic spin fluctuations. Near cubic sym- 
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FIG. 1: Numerical systomatics in the evaluation of pair- 
ing amplitude, rup, versus the scaled temperature, fcsT/Wo, 
where Wo is the non-interacting electron bandwidth, for two- 
dimensional lattices with 0.85 electrons/site and U/Wo = 0.5. 
Top: Variation as a function of number of Matsubara frequen- 
cies, m. Middle: Dependence on the number of lattice points, 
Nl- Bottom graph: Variation as a function of DCA cluster 
size, Nc. 
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FIG. 2: Fluctuation exchange approximation results for the 
mean-field superconducting phase diagram of the Hubbard 
model as a function of the ratio of the interplanar to intrapla- 

nar hopping, tz/txy, and reduced temperature, ksT/Wo, for 
0.85 electrons/site and U/Wo = 0.5. The d.^2_y2 and d^^2_r2 
pairing states arc most stable for tz/txy < 1 and tz/txy > 1 
respectively. For tz/txy ~ 1 a second transition occurs as 
a function of decreasing temperature into the time-reversal 
symmetry breaking state d^2_y2 + id^^2_^2. 



metry, two phase transitions are realized, the first corre- 
sponds to a transition to a single-component state and 
the second transition occurs to a multicomponent state 
with broken time reversal symmetry. On the basis of 



these calculations, we propose that Cclna is a candidate 
material for exhibiting broken time-reversal symmetry in 
the superconducting state. 
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